Nonlocal transport through multiterminal diffusive superconducting nanostructures 
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6J[). Motivated by recent experiments on nonlocal transport through multiterminal superconducting 

hybrid structures, we present self-consistent calculations based on quasiclassical Green's functions 
for the order parameter, currents and voltages in a system consisting of a diffusive superconductor 
connected to two normal and one superconducting electrodes. We investigate non-equilibrium effects 
for different biasing conditions corresponding to measurements of the nonlocal conductance and of 
the nonlocal resistance. It is shown that while the nonlocal conductance does not change its sign, this 
change might be observed in a nonlocal resistance measurement for certain parameter range. The 
change of sign of the nonlocal signal takes places at a voltage of the order of the self-consistent gap 
' of the superconducting region. We show that this is not related to the nonlocal Andreev processes 

but rather to non-equilibrium effects. We finally discuss the case of four terminal measurements and 
demonstrate that a change of sign in the nonlocal resistance appears when the current injected into 
the superconductor exceeds a critical value. The connection to the existing experiments is discussed. 



I. INTRODUCTION 

The possibility to create and control entangled electron pairs from a superconductor has renewed the interest in 
transport through superconductor-normal metal hybrid structure a 1 ' 2 ' 3 ' 4 ' 5 . The basic idea is to exploit the long-range 
coherence of the Andreev reflection^, in order to couple two spatially separated normal electrodes connected to a 
O . superconducting region within a distance of the order of the superconducting coherence length £5. An Andreev 
processes that takes place at two different interfaces is called a "crossed"' Andreev reflection (CAR). As suggested 
in Refp£ spatial correlations can be probed by nonlocal transport experiments. A typical setup for detection of CAR 
processes consists of a grounded superconducting region (S) connected to normal electrodes (N). The information on 
such processes would be encoded, for instance, in the voltage which is measured in one of the S/N interfaces when 
a current is injected through the other one. Besides CAR processes, individual electrons can also tunnel across the 
superconductor. This normal tunneling has been called "elastic cotunneling" (EC). CAR and EC contributions to the 
non-local conductance have opposite signs, and in the lowest order of tunneling cancel each other—. For higher orders 
in the tunneling the EC dominates over CAR 9 - and the nonlocal signal becomes finite. Surprisingly a change of sign in 
the nonlocal resistance and conductance was reported as a function of the local voltage for a NSN layered structure 2 - 
and a FSF multiterminal structure 4 (F denotes a ferromagnetic metal). According to Refi 2 - by low (high) voltages EC 
(CAR) processes dominate the nonlocal transport. Similar behavior was observed in Ref4. for samples with high S/F 
barrier resistance. The latter experiment also showed dominance of CAR processes at low voltages for samples with 
higher interface transparency. These experimental results have lead to several theoretical work a 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 , 
which attempt to find a microscopic description for those observations. However, up to now theories based on non- 
interacting models could not explain the change of sign of the non-local conductanc e 10 ' 11 ! 12 ' 13 . For a layered NSN 
structure, as the one of Refi 2 -, the change of sign of the non-local conductance has been explained by taking into 
account interaction of the conducting electrons with their electromagnetic environment 1 ^. The later description is 
valid in the tunneling limit, thus the observation of negative non-local conductances in the case of good interface 
transparencies remains without microscopic explanation yet. 

In principle, a description of the non-local transport in terms of CAR and EC is only valid in the tunneling limit. 
In this case the system is in a quasi-equilibrium state, i.e. the current (or the corresponding bias voltage) is much 
smaller than its critical value and the distribution function of quasiparticles is the equilibrium one. This assumption 
has been made in most of the theoretical works mentioned above. In particular the superconducting gap was assumed 
to have the bulk value, i.e. the superconducting order parameter was not affected by the possible deviation of the 
distribution function from its equilibrium value. In some experiments though, this is not the case. For example in 
Refsi 1 -' 3 -^ the transparencies of the S/N interfaces are not necessarily low, and in Refi 3 - the current injected into the 
S region reached its critical value. Thus, for a proper description of these experiments one needs to go beyond the 
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FIG. 1: (Color online.) Sketch of a three-terminal device for (a) measurement of the non-local conductance and (b) measurement 
of the non-local resistance, (c) Typical experimental setup for measurement of the non-local resistance. 

quasi-equilibrium approach. A first attempt was done in Refai^. It was shown that non-equilibrium effects may 
play a crucial role on the non-local transport properties. In particular, in Ref self-consistent calculations based on 
a two dimensional tight binding model were implemented. It was shown that far from the quasiequilibrium regime the 
nonlocal transport cannot longer be described in simple terms of EC and CAR processes. For some set of parameters, 
a change of sign in the non-local resistance was obtained. This change of sign is not related to the predominance of 
CAR but rather to the possibility of having a negative local conductance at the interface where the current is injected. 
In RefJ^ a non-monotonic behavior of the non-local resistance as a function of the temperature was obtained, which 
resembles the observations of RefA However, the non-local resistance as a function of the injected current or bias 
voltage was not investigated in that work. 

In the above mentioned theoretical works, one computes the non-local conductance in a three terminal device. In 
other words, one assumes that one of the normal terminals is biased to a voltage Vl, while the second normal terminal 
is grounded. In this way one determines the current Ir flowing into the latter terminal and computes the non-local 
conductance G n i = cHr/cIVl. Experimentally, however, it is simpler to fix the injected current through one of the 
normal terminals and measure the voltage induced at the second normal terminal where no current is flowing. Thus, 
the measured quantity is the non-local resistance. Moreover, some experiments were performed in a multi-terminal 
geometry^. 

In this paper we present a complete self-consistent theory for the nonlocal transport through a diffusive supercon- 
ducting region connected to several normal electrodes. We calculate both, the non-local conductance and resistance 
with the help of the quasiclassical Green's functions (GF) approach. In a first part we concentrate on a three terminal 
device, where the current is injected from a normal electrode Nl into the superconducting region S, maintaining 
the second normal electrode Nr grounded. We determine the self-consistent gap, the current flowing into Nr and 
compute the non-local conductance G n i. We show that our model, as in previous works, predicts no change of sign 
for G n i. In a second part we consider again a three terminal device, but now we assume that no current is flowing at 
the S/Nr interface terminal. We then determine the self-consistent gap, the voltage induced in Nr, and determine 
the non-local resistance R n i . For certain range of parameters we obtain a change of sign R n i due to the appearance of 
a negative local conductance at the N^/S interface. Finally, we considered a four terminal setup. Again we assume 
that the current through the S/Nr interface is zero, but now we determine the voltage induced in the Nr electrode 
measured with respect to the end of the superconductor in which no current is flowing. In this case we obtain a change 
of sign as in the experiment of RefA As we show below the origin of this change of sign is not due to a negative local 
conductance, but to the non-equilibrium distribution created in the superconductor by the injected current, which 
eventually leads to a transition into the normal state. 

The rest of the paper is organized as follows: in section II we introduce the model and basic equations used 
throughout the paper, section III and IV are devoted to analyze the results for the three terminal situation (we 
discuss the conductance measurement conditions in section III and the resistance measurement case in section IV), 
and finally in section V we analyze the four terminal case. Some concluding remarks are given in section VI. 
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II. THE MODEL AND BASIC EQUATIONS 



A typical experimental setup for the measurement of non-local transport properties is shown in Fig. lc. On top 
of a nanoscale superconducting wire one places several (in our case two) normal wires. A current is injected from 
one of the normal electrodes and flows in the direction shown in the figure. The non-local resistance is then obtained 
by measuring the potential difference between the other normal electrode and the end of the superconducting wire 
through which no current is flowing. A strong enough non-local signal is measured when the distance between the 
normal leads does not exceed much the coherence length of the superconductor. 

A self-consistent computation of the spatial variation of the gap and currents in such a device is a formidable task 
which we will not address here. Nevertheless, the mean features of the system can be described by considering the 
geometry shown in Figs, la and lb. A superconducting region S is connected to two normal electrodes Nl and Nr 
and to a superconducting electrode S e , made of the same superconducting material as S. The three interfaces will 
be described by characteristic energies th^R and es, defined below. As we are not interested in the spatial variation 
of the non-local correlations, we simplify the problem by assuming that the central superconducting region S has 
dimensions smaller than the superconducting coherence length, which in the diffusive limit is given by £g = y/D/A.. 
Here A is the superconducting gap and D is the diffusion coefficient. Thus, we may assume that the order parameter 
and the non-equilibrium distribution are uniform over 5™". 

In order to determine either the non-local conductance (Fig. la) or resistance (Fig. lb) we need to calculate 
the current density j and the self-consistent order parameter A in the S region. As one can see from our results 
below, we always find a stationary current state. This behavior is distinctive for junctions consisting of a mesoscopic 
superconductor in contact with a bulk one2£. Thus, j and A can be expressed in terms of the Keldysh component of 
the matrix GF g 

j = ^/*Tr{^}-, (2) 

where A is the BCS coupling constant which determines the critical temperature, and Rn is the normal state resistance 
of the S region. The function g is a 4x4 matrix in the Nambu®Keldysh space with the usual structure 



g R g K 

g A 



(3) 



while g are 2x2 matrices in Nambu space. In the diffusive limit these functions are the solutions of the Usadel 



equation^ 



DV(gVg) - ie [f 3 , g] - i [A, g] = -i [t in ,g] , (4) 



supplemented by the normalization condition g 2 — 1. Here E,„ is the self-energy term describing inelastic processes. 
In the time relaxation approach Sj n is proportional to 1/Tj n , where Ti n is the inelastic relaxation time. We will 
assume that 1 /r in is the smallest energy scale and neglect this term. At the interfaces with the electrodes we use the 
Kupryianov-Lukichev boundary conditions^ 

DgVg\ n = eid[g,gi] , (5) 

where & are the GF of the electrodes(i = L, R, S e ), ej = erh/^rB, erh = D/d 2 is the Thouless energy, tb = Rbi/Rn, 
Rbi is the ith barrier resistance per unit area, and n denotes a unit vector normal to the interface. We assume 
that the GFs of the electrodes remain unchanged and equal to the bulk values, i.e. g R ^ = ±t 3 in the normal 
leads, and g^ A ^ = g'scs^ + Ibcs^ 2 m ^ ne superconductor electrode, where g^^g — e/y/(e± irf) 2 — Aq and 
f]3cs — ^-o/ i ir/) 2 — Aq. While the Keldysh components are given by 



- g*Pi F t g A (6) 



where 



F l+ T + Fi^T 3 , 



Vi is the voltage in electrode i , and Fi± — | [tanh( e tS?* ) ± tajih( e "ZS?* )] . We also assume that Vs e = 0. In principle 
the boundary conditions Eq. ^ are valid for low transmitting interfaces. In the present work we consider that the 
interface transparencies may vary in the range 10 -3 — 10 _1 for which Eq. ([5]) is sufficiently reliable. 
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With the help of Eq. ([5]) we can calculate the total current at each interface using the expression: 

eIiR N = J deTrr 3 [g i: g] K . (7) 

In this case the GF inside S does not vary considerably and the Usadel equation (|3]) can be integrated over space 
coordinates using the boundary conditions Eq. (0 . In this way one obtains a set of algebraic equations which can be 
written in a compact form 

[A,s]=0 (8) 

where 

A = ^ e i9i + er 3 + A - E in . 

i=L,R,S 

Eq. ([5]) is equivalent to the Nazarov's circuit theory equations 22 , which were used in Refs i n i 23 for nonlocal transport 
calculations. The solution for the R,A and K components of g which satisfy Eq. ([8]) and the normalization condition 
can be formally be written as: 



*R{A) 



(9) 



-k _ k K -g R A K g 



Va r a r + Va a a a 



(10) 



Substituting these expressions into Eqs. (jT][3]) enable us to obtain numerically the self-consistent order parameter, the 
currents through the interfaces and the non-local voltage induced at the right electrode in the resistance measurement 




FIG. 2: (Color online). The amplitude of the self-consistent order parameter, the current II injected from the left normal 
electrode and the current Ir measured at the right electrode, as a function of the voltage Vl- Panels (a)-(c) for es ~ 0.2, 
el = 0.1, T — O.OIAq and different values of en. Panels (d)-(f) for es = 0.2, e_R = 0.1, T — O.OlAo and different values of cl- 
We have defined R = RxeTh/Ao- 



III. MEASUREMENT OF THE NON-LOCAL CONDUCTANCE 



In this section we consider the experimental set-up of Fig. la. The left normal electrode is biased by a voltage source 
at Vl. For a non-local conductance measurement we will assume that the Nr is grounded (Vr = 0), and compute the 
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FIG. 3: (Color online). The amplitude of the self-consistent order parameter as a function of the voltage Vl, for different 
temperatures. Notice that in the range of temperatures shown and for large values of Vl the gap is enhanced by increasing T. 



current Ir through the interface S/Nr from Eqs. (|7l9ll0p and the self-consistent order parameter A from Eq. (TJ). 
Before we address the non-local properties of the system let us discuss the results concerning local properties. In Fig. 

we show the amplitude of A as a function of the bias voltage Vl for es = 0.2Ao, (l = O.lArj and three different 
values of €r = 0.01 — 0.1, corresponding to transmission coefficients in the range 10~ 3 — 10~ 2 (we assume that the 
length of the superconducting region is about bOnm). All energies are given in units of Ao, which is the value of the 
order parameter in the bulk at T = 0. For a fixed low value of Vl, A is reduced by increasing the coupling with 
the right normal electrode, which is a consequence of the inverse proximity effect. At some value V£ ~ O.8A0 of the 
order of the self-consistent A, one can see an abrupt reduction of the self-consistent order parameter. For voltages 
larger than V£ the quasi-particle current through S becomes considerably larger (Figs. [2ja and[2j;), i.e. the system 
is driven out of equilibrium. An interesting consequence of this non-equilibrium state for voltages Vl >~ V£ is the 
enhancement of the self-consistent gap by increasing the temperature, as shown in Fig. [3] This effect is related to 
the stimulation of superconductivity by quasiparticle currents in SIS systems, and was studied both theoretical^ and 
experimentally^. 

From Fig. [^1 one can also see that the suppression of A at V£ becomes more abrupt the weaker the coupling with 
the right electrode is. Figures [2Jd and also show the corresponding current at the left and right interfaces. As 
expected the larger the resistance of the right interface (small cr) the smaller the value of Jr. A strong non-equilibrium 
situation takes place when most of the current injected flows into the superconducting electrode, i.e. when €r is small 
enough (in our example = 0.01). In this case the gap becomes multivalued and this is reflected in the behavior 
of the currents II and Ir. Multivalued solutions for the self-consistent gap were also found in Ref*2!i for SIS systems 
and recently in Refill for a NSN system. 

Another interesting feature of this system is the existence of a region of voltages for which the local conductance 
is negative (see Fig. [5Vj). This behavior was also obtained in Ref.~ for a system consisting of a superconducting link 
separating a normal and a superconducting electrode. Also in Ref. 14 negative local conductance was obtained for a 
two dimensional ballistic superconductor attached to two normal electrodes. Notice however, that the part of the 
curve corresponding to a negative conductance would be not accessible in current biased experiments. If we now fix 
the value of 6r at its maximum value cr = 0.1 and vary €l, we see that even for the smallest coupling (ej, = O.OlAo) 
the variation of the gap and the currents is smooth and no signatures of multivalued solutions appears for this range 
of parameters (see bottom row of Fig. [2]) . 

Let us now focus on the non-local transport and compute the non-local conductance, which can be obtained easily 
from the knowledge of Ir (Fig. [^fc and [2f ) . It is given by the expression 

ffi, - gf ■ (11) 

and shown in Fig. 0] as a function of the bias voltage Vl, for different values of the coupling energies cl r- In all cases 
G n i is very small in the region of low voltages. This is in agreement with the zero non-local conductance obtained 
in the lowest order of tunneling^, due to the cancellation of the EC and CAR processes. However, our results are 
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FIG. 4: (Color online). The voltage dependence of the (negative) non-local conductance G„i normalized with respect to 
R = R N e T h/A . We have chosen T = 0.01A , e s = 0.2, e L = 0.1, T = 0.01A and different values of e R . 



in all order of tunneling and therefore a dominance of EC (negative G n i) is obtained in accordance to RefA For 
voltages of the order of the self-consistent A, the value of G n i becomes significant. However, no change of sign is 
observed. As mentioned in the introduction, non-interacting models do not exhibit a change of sign of the non-local 
conductance in three-terminal NSN structures. This can only be achieved when electron-electron interactions are 
taken into account^. 



IV. MEASUREMENT OF THE NON-LOCAL RESISTANCE IN A THREE-TERMINAL DEVICE 

In a real experiment it is easier to measure a voltage rather than a current. Indeed the experiments of Refsi 2 -'^'^ 
were performed in (a) the current biased regime and (b) instead of the current, the non-local resistance (or voltage) 
was measured. Theoretically, it is not simple to impose a current bias. Therefore we will still work in the voltage 
biased case but determine the induced non-local voltage and resistance imposing zero current at the S/Nr interface. 
Thus, all the current injected from the left normal electrode flows into the superconducting electrode S e (see Fig. lb). 
The current at the right interface is given by (c/. Eq.([7])): 

eI R R BR = \J deTrg K (e) - ~ J dev a {e) (tanh(i±|^) - tanh(^|^)) , (12) 

where v s (e) — (g R — g A )/2 is the density of states of S and Rbr is the R barrier resistance per unit area. The first 
term in the r.h.s is proportional to the quantity Q* identified in the literature as the charge imbalance potential 27 , 
which appears due to a non-equilibrium distribution in the superconductor. The second term is the usual quasiparticle 
current term. The voltage Vr is measured with respect to the ground (see Fig. lb) and it is obtained by imposing 
Ir = Q. 

The results for A, the injected current II and the induced voltage Vr as a function of the bias voltage Vl are shown 
in Fig. O for fixed values of es,L and different values of €r in the same range as in Fig. [51 In the case of low values 
of cr, the self-consistent gap has a very similar behavior as in the preceding section. However, for the largest value 
€r = O.lAo the suppression of A is larger as the one obtained by imposing Vr = 0. Notice also that the region of 
negative local conductance (Fig. [SJd), associated with the abrupt change of the order parameter appears now for all 
values of er. It is clear that by having imposed Ir = all the current injected must flow through the S/S e interface 
and our system behaves similarly to the N/S/ S studied in Ref*2&. In Fig. [SJ: we also show the voltage Vr induced in 
the right electrode calculated by equalizing (fT2"| to zero. At low Vl values the induced voltage Vr is very small, but 
it experiences a jump at the value of Vl where the self-consistent gap exhibits its maximal drop. 

We have now all quantities required to compute the non-local resistance which would be measured. This can be 
calculated from the expression: 



= dVx = dVx (dI L \ 
1 dI L dV L \dV L J 



(13) 
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The measured voltage Vx depends on the experimental set-up. We are considering here the three-terminal structure of 
Fig. lb, and determining Vji respect to the ground. Thus in this case Vx = Vr, i.e. the one shown in\5jp. According to 
Eq. (|13p there are two factors determining the non-local resistance. One which is the inverse of the local conductance 
^ll = dVL,/dlL and which decrease by increasing €l- The second contribution is given by dVn/dVL which is nonzero 
only if a non-equilibrium distribution appears in the S region and is related to the charge imbalance term (c/. Eq. 

«■ 

In Fig. [5] we show the dependence RniiYh) for two values of €r and = 0.1A , £s = 0.2A , T = 0.01A . The 
change of sign of R n i is a consequence of the negative local conductance which appears between the values Vl — O.8A0 
and Vl — Ao (see Fig. [5b). The position of the first peak of Rni is determined by the value of Vl at which A drops 
substantially, while the position of second peak is determine by Ao. As mentioned above, in current bias experiments 
the curve R n i may look very different to the ones shown in Fig. [6l since the region of negative local conductance may 
not be observed. 

In experiments as those of Refsi^^ the current flows along a superconducting wire, while in our model it flows into 
the reservoir S e . For simplicity we have assume that the latter remains unaltered for all values Vl considered here. 
In particular the value of the gap is the bulk BCS one for any value of Vl- However, in the experiments, when the 
current flowing through the wire reaches the critical value the superconducting gap is suppressed homogeneously in 
the region where the current is flowing. This leads to the observation of only one peak in the non-local resistance^. 
In the next section we will model this situation. We should also emphasize that the change of sign of the non-local 
resistance obtained in Fig. [6] is due to the fact that the local conductance Gll = AIlIAVl becomes negative for some 
values of Vl (cf. Fig. [5b). If the coupling es is large enough, the local conductance remains always positive and so the 
non-local resistance. This is shown in Fig[7| where the amplitude of the self-consistency gap, the current injected and 
the voltage induced at the right electrode are plotted as a function of Vl, for 65 = 1, 1.5 and 2Ao- If one compares 
these results with those obtained for a smaller £5 (Fig. [5J, one sees that A is now only weakly suppressed and that 
the current II increases monotonically. Thus the non-local resistance is always positive as it shown in Fig. [5J 

Finally, we show in Fig. [5] the temperature dependence of the zero bias non-local resistance for different values 
of the coupling parameter £5. We can see that while for the small values of £5 the non-local resistance increases 
monotonously with the temperature, for larger values of es, R n i reaches a maximum value. In the latter case the 
charge imbalance effect becomes important and dominates over the local conductance factor for large temperatures. 
This behavior is in agreement with previous calculations of R n i in a superconducting quantum dot^. One could 
conclude as in Re£i^ that the non-monotonic behavior of R n i (T) is in qualitative agreement with the observations of 
Refsfi^. However, we hardly believe that. If this would be the case then one should obtain for the same parameter 
range a change of sign for R n i as a function of the injected current, as observed in the experiments^. On the contrary, 
Fig. [5] clearly shows a monotonic increase of R n i with the applied voltage. In the next section we will show that the 
change of sign of i? n ; observed in the experiment is due to the suppression of the superconductivity by the injection 
of a current. Also the peak of R n i observed as a function of temperature could be understood within this model. 
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FIG. 6: (Color online). The voltage dependence of the non-local resistance calculated from Eq. (|13p for es = 0.2Ao, el = O.lAo, 
T = O.OIAq, and two values of ej{. 



V. FOUR TERMINAL STRUCTURE 



We now consider a situation closer to that of the experiments of Reb^i, in which the non-local resistance (volt- 
age) has been measured in a multi-terminal setup consisting of a superconducting wire attached to several normal 
terminals. In these experiments the nonlocal voltage corresponds to the potential difference between one end of the 
superconducting wire and one of the normal leads. We model these experimental situations as shown in Fig. I10I On 
the top of a superconducting wire we place two normal contacts. A current flows from the Nl contact to the N e due 
to the bias voltage Vl applied between the contacts. We are interested in the voltage difference measured between 
the end of the wire, which we denote by S', and an additional normal contact Nr. We use on purpose the same 
notation as in Fig. 1 in order to use straightforwardly the expressions derived in section II. The only difference is that 
the drain electrode is now in the normal state and that we measure the voltage difference between the end S' of the 
superconducting slab and the normal electrode Nr. The electrode N e is grounded. In order to compute the resistance 
measured between S' and Nr we proceed as in the last section, determining the Green functions from Eqs. (j51[TT)|) 
and the self-consistent gap. Now we impose that the currents through the S/Nr and through the S/S' interfaces are 
zero. The current trough the S/ S' interface can be written as the sum 

I S /S' = I J + Iqp (14) 

of the Josephson and the quasiparticlc contribution respectively The first is given by the product of anomalous 
Green's functions in Eq. ^ while the latter by the product of normal components. Since no current is flowing into 
S' we assume that there the Green's functions are those in equilibrium, with |A| equal to the self-consistent |A| in 
the S region at Vl = 0, and a phase <fi which is determined by imposing the condition of no current through the S/S' 
interface. For voltages lower than a critical voltage V£ ~ A we always found a finite value of (j>. In this case the 
voltage V induced in S' equals zero, and the measured V n i coincides with Vr of the previous section. The voltage 
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eV,/A n eVJ\ eV,/A n 

FIG. 7: (Color online). The amplitude of the self-consistent order parameter, the current II injected from the left normal 
electrode and the induced voltage Vr measured at the right electrode, as a function of the voltage Vl, for e_L = = 0.1, 
T — O.OIAq and different values of es- 




eV L /A Q 

FIG. 8: (Color online). The voltage dependence of non-local resistance for the same parameters as in Fig. 0. 

V£ is the voltage at which the self-consistent gap vanishes, i.e. when the current flowing in S reaches its critical 
value. For values of Vl larger than this value the quasiparticle current becomes finite, and a voltage V is induced 
in 5". We compute it from an expression obtained by equalizing Eq.([7]) to zero. Thus, the non-local resistance is 
given by Eq. ([13]) where now V x = V R for V L < VI and V x = Vr - V for Vl > V£. In Fig. [TT] we show the 
result of our calculation for R n i as a function of Vl- It shows a peak at V£. Since the latter is of the order of A, the 
peak is shifted to lower voltages by increasing the temperature (Fig. ITT]) . This behavior is in agreement with the 
experimental observation of Ref£, where the peak occurred at values of the bias current close to the value of critical 
current of the superconducting wire. The change of sign of R n i is related to the non-equilibrium situation created in 
S by the injection of a current from Nl- We also show in the inset of Fig. [TTJthe temperature dependence of R n i, 
which exhibits a pronounced peak for Vl — V£. This is again in agreement with the observations of Ref4. These 
results demonstrate that the model studied here contains the main ingredients for describing experiments on non-local 
transport as the one of RefA Within this model the change of sign of the non-local resistance has its origin in the 
deviation of the distribution function of the superconductor from the equilibrium one. Notice, that as in Rcf. 3 , the 
change of sign of R n i occurs at the critical current which corresponds to the voltage V£ in our model. 

In other experiments^, however, the change of sign occurred at lower voltages. This discrepancy is at the moment 
not clarified and may be related to the inclusion of electron-electron interactions as proposed in RefJ^. 
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FIG. 10: (Color online). Sketch of the 4-terminal structure under consideration. 



VI. CONCLUSIONS 



We have presented a self-consistent analysis of the transport properties of a structure consisting of a mesoscopic 
superconductor whose dimensions are smaller than the characteristic length £g, attached to two normal and one 
superconducting terminals. We have analyzed two measurement methods: one in which the detector (R) electrode is 
grounded and the leaking current is measured, and one in which the current through this lead is fixed to zero and 
the induced voltage is measured. In both cases we observe that the self-consistent order parameter in the mesoscopic 
central region exhibits an abrupt drop at a certain voltage of the order of the self-consistent A. Associated to this 
drop the local differential conductance at the injector lead (L) may become negative for certain values of the coupling 
parameters, resulting in a change of sign of the non-local resistance. As we stress throughout this manuscript, this 
change of sign would not be related to a dominance of CAR over EC processes but to a non-equilibrium effect. We have 
still described another mechanism for the appearance of negative non-local resistance which is probably most suitable 
for explaining the observations of Ref. 3 . This mechanism is applicable in a four terminal geometry and corresponds to 
the injection of large currents which may switch the superconducting region into the normal state. The observation of 
a change of sign in the non-local signal at smaller bias and small transparencies like the ones reported in Refsi^£ are 
certainly not possible to be explained with the theoretical model presented in this work and might be related to the 
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FIG. 11: (Color online). The non-local resistance measured in a four terminal structure as a function of the bias voltage Vt 
for e_L = e_R = O.lAo, es = 0.2, and different values of the temperature. Inset: The non-local resistance as a function of the 
temperature at zero bias. 

influence of electron-electron interactions, as already pointed out in Ref.— . Further work for analyzing the combined 
effect of interactions and non-equilibrium effects is under progress. 
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